Persistent currents in Bose gases confined in annular traps 
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We examine the problem of stability of persistent currents in a mixture of two Bose gases trapped 
in an annular potential. We evaluate the critical coupling for metastability in the transition from 
quasi-one to two-dimensional motion. We also evaluate the critical coupling for metastability in 
a mixture of two species as function of the population imbalance. The stability of the currents is 
shown to be sensitive to the deviation from one-dimensional motion. 
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I. INTRODUCTION 

Bose-Einstein condensates of dilute vapors of atoms 
offer a very promising testing ground for questions as- 
sociated with superfluidity for a number of reasons. 
Firstly, these gases are dilute as opposed, for example, to 
the "traditional" superfluid liquid helium. Furthermore, 
the atomic gases can be manipulated in many different 
ways, including the shape of the confining potential, the 
strength and the sign of the effective interatomic interac- 
tion, the number of different species in multicomponent 
systems, etc. 

While the term "superfluidity" covers a whole collec- 
tion of many different phenomena (l|, we focus in the 
present study on the metastability of superflow in an- 
nular traps 043- Persistent flow has been observed re- 
cently in a Bose-Einstein condensate of sodium atoms 
confined in a toroidal trap [J]. In this experiment, an 
initial angular momentum of % per particle was trans- 
ferred to the atoms and the rotational flow was observed 
to persist for up to ten seconds, limited only by the trap 
lifetime and other experimental imperfections. Persistent 
currents with two units of angular momentum were also 
observed in the same experiment. 

Theoretical studies have examined the existence of 
metastable states in Bose-Einstein condensed gases 
trapped in single- 04§| and double-ring-like confining 
potentials, as well as other phenomena in double-ring-like 
traps [111, EH • As shown in Ref . [|| , mixtures of two com- 
ponents that are confined in a strictly one-dimensional 
single ring support persistent currents, but the interac- 
tion strength necessary for metastability increases with 
the admixture of the second component. For compara- 
ble populations of the two components, the critical value 
of the coupling becomes infinite. Finally, persistent cur- 
rents with a value of the circulation higher than one unit 
were shown to be stable only in single-component sys- 
tems. 

In the present study we investigate the effect of the 
deviations from purely one-dimensional motion on the 
metastability of the currents. We also examine two- 
component Bose-Einstein condensates, and find that it 
is possible to have metastable superflow for sufficiently 
small admixtures of the second component. One of the 



novel results of our study is that in this latter case the 
deviation from one-dimensional motion gives rise to per- 
sistent currents with circulation higher than one unit, as 
opposed to the purely one-dimensional case. 

The paper is organized as follows. We describe in 
Sec. II our model, and in Sec. Ill the two methods that 
we have employed to solve it, i.e., the mean-field approx- 
imation and numerical diagonalization of the Hamilto- 
nian. In Sec. IV we present our results. More specifically, 
in Sec. IV. A we consider the single-component case and 
study the existence of metastable states in the transition 
from quasi-one-dimensional to two-dimensional motion. 
In Sec. IV. B we examine the effect of the admixture of a 
second component on the stability of the persistent cur- 
rents, comparing the obtained results with the ones of 
strictly one-dimensional motion. Finally, we present our 
summary and conclusions in Sec. V. 



II. MODEL 

Let us consider two distinguishable species of bosonic 
atoms, labelled as A and B, with populations Na and 
Nb respectively; without loss of generality we assume 
below that Nb < Na- We also assume two-dimensional 
motion, which in an experiment corresponds to the case 
of very tight confinement in the perpendicular direction. 
In the plane of motion we model the annular trap by a 
displaced harmonic potential, 



V{p) = \muj 2 { P - R ) 2 , 



(1) 



which is plotted in Fig. [TJ Here p is the usual radial 
variable in cylindrical polar coordinates, i?o is the radius 
of the annulus, ui is the trap frequency, and M is the 
atom mass, which we assume for simplicity to be equal 
for both species. We stress that the results presented 
below with the potential of Eq. (JTJ are - at least qualita- 
tively - similar to the ones that we have obtained using 
a harmonic-plus-Gaussian trapping potential. 

In order to investigate the effects due to the transition 
from quasi-one-dimensional to two-dimensional motion, 
we vary the confinement strength ui in the radial direc- 
tion. For large w, the confinement is strong enough to 
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eigenenergies and the corresponding eigenvectors. This 
can be done for a range of the values of L, yielding the 
dispersion relation E(L), i.e., the smallest eigenvalue for 
some given L, allowing us to identify the possible local 
minima associated with the metastable states that give 
rise to persistent currents in the system. 

The basis states that we choose in this approach are 
the eigenfunctions of the single-particle problem, which 
are of the form 



(4) 



FIG. 1: (Color online) The confining potential V of Eq. (0 as 
a function of the cartesian coordinates (x,y), for Ro/ao = 4 
and lj/ujq = 1. 



freeze out the motion in this direction, which makes the 
system quasi-one-dimensional. As u> decreases, the con- 
finement becomes weaker and, as a result, the motion be- 
comes two-dimensional. In all the calculations presented 
below we fix the radius of the annulus to Ro/ao = 4, 
where a — (Moj a ) is the ususal oscillator length 
corresponding to some fixed frequency lu . 

The interatomic interactions are modelled via the usual 
effective contact potential, which is assumed to be re- 
pulsive. Therefore, the Hamiltonian H of the system is 
H = H sp + Hi n t , where H sp is the single-particle part 
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and -ffi n t is the interaction part, given by 
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N A ,N B 
t=l,J=l 



(3) 



Here the parameters Uki are proportional to the s-wave 
scattering lengths for zero-energy elastic atom-atom col- 
lisions. For simplicity, we take uaa = ubb = uab = u, 
which is also experimentally relevant in several cases; 
furthermore, since the interaction is assumed repulsive, 
u > 0. 



III. METHOD 

We attack this problem using both the mean-field 
Gross-Pitaevskii approximation, as well as diagonaliza- 
tion of the many-body Hamiltonian. 

Within the diagonalization approach, we fix the parti- 
cle numbers Na and Nb and the total angular momen- 
tum Lh, and use the Lanczos method 11311 to obtain the 



The quantum number n corresponds to the radial excita- 
tions, m is the quantum number associated with the an- 
gular momentum, which results from the rotational sym- 
metry of our problem, and 9 is the usual angular variable 
in cylindrical polar coordinates. The radial wavefunc- 
tions R n ,m(p) are evaluated numerically. If E n ^ rn are the 
corresponding eigenvalues, we work under the assump- 
tion that the typical value of the interaction energy Vint 
is weak enough, so that V^t "C Si,o- Thus, we set n = 
and the only quantum number that remains is m, whose 
highest value is chosen so that V ln t ^ £t),m "C £a,o- Hav- 
ing evaluated the single-particle states, these are then 
combined in all possible ways to form a basis of the Fock 
states for the many-body problem for some fixed Na, 
Nb, and L. 

Turning to the mean-field approximation, within this 
approach the system is described via the two order pa- 
rameters of the two components, <Pa and <f>s- From 
Eqs. ([2]) and ((3]) one obtains the two coupled Gross- 
Pitaevskii-like equations for these two order parameters 
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*>a + V{r)cf>A + 9aa\(I>a\ 2 (I>a + 9ab\4>b\ 2 (I)a 



>B + V{r)(f>B + gBB\(t)B\ 2 (f>B + 9Ab\4>a\ 2 (I)B 



Ha4>a, 

ia\ 2 4>b 
Vb4>b- (5) 



In the above equations the two order parameters are nor- 
malized as /|0A| 2 GPr = 1, and J \(j) B \ 2 d 2 r = N B /N A - 
The parameters gij are defined as qaa = Nauaa, 9bb = 
Naubb, and qab — Nauab- Finally, ha and \xb are the 
chemical potentials of the two components. 

To solve Eqs. ([5]) we make use of a fourth-order split- 
step Fourier method within an imaginary-time propa- 
gation approach Starting with an initial state 
with some given angular momentum, the imaginary-time 
propagation drives the system along the dispersion rela- 
tion E(L) until it finds a local minimum. If E(L) has such 
local minima, the final state may be a metastable one. 
Otherwise, the initial state decays to the non-rotating 
ground state of the system. One should note that the 
dissipation of the angular momentum L is inherent to the 
diffusive character of this method, and it does not con- 
tradict the fact that the Hamiltonian commutes with the 
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operator of the angular momentum, since the imaginary- 
time propagation operator is not unitary, and therefore 
does not satisfy Heisenberg's equation of motion [jjj . 

The two methods mentioned above complement each 
other: the diagonalization gives the full dispersion rela- 
tion E(L) for any fixed angular momentum Lh, but due 
to numerical limitations in diagonalizing matrices of large 
size, one has to restrict the number of particles, as well 
as the values of the interaction strength. The mean-field 
Gross-Pitaevskii approach, on the other hand, does not 
have these limitations, but does not allow us to evaluate 
the whole dispersion relation E(L) - at least in a straight- 
forward way, but rather its local/absolute minima. Fi- 
nally, due to its mean-field character, possible correla- 
tions between the atoms in certain limiting cases are not 
captured within this approximation, although they are 
not expected to be important in the problem that we 
consider here. 



IV. RESULTS 

A. Deviation from one-dimensional motion and 
metastability 

We start with the case of a single component, i.e., when 
the total number of atoms ./V = Na + Nb is equal to Na- 
We evaluate the minimum interaction strength g mm nec- 
essary for the system to support persistent currents as 
a function of the width of the annulus, or equivalently 
as a function of the confinement frequency uj. As we ar- 
gued earlier, for large values of uj, the motion becomes 
quasi-one-dimensional, while as ui decreases, the motion 
becomes two-dimensional. We solve Eqs. (J5J) and also 
diagonalize the Hamiltonian numerically for several val- 
ues of uj. The results of both calculations are shown in 
Fig. [21 represented by (black) circles (mean-field) and 
(red) crosses (diagonalization). As seen from this plot, 
the results obtained from the two approaches are in good 
agreement with each other. 

At this point it is instructive to see how these data 
compare with the analytical result of a purely one- 
dimensional model. To reduce the two-dimensional prob- 
lem into an effectively one-dimensional one, at least in 
the limit where the width of the annulus (set by the 
oscillator length) is much smaller than its radius Ro, 
one may start from the initial expression for the en- 
ergy and integrate over the radial degrees of freedom 
|16| . To simplify this calculation, we assume that the 
density in the transverse direction is nonzero only for 
Ro — flosc/2 < p < Rq + a osc /2, where Rq is the radius of 
the annulus and a osc is the oscillator length correspond- 
ing to uj. In this limit, the order parameter may be writ- 
ten in a product form, 



dimensional nonlinear equation 



4>a(p,6) = Ra(p) &a(0). 



(6) 
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where the coefficient of the nonlinear term involves the 
inte gra l J \Ra(p)\ 4 pdp. According to previous studies 
[l7l Il8| , Eq. ([7]) implies that the minimum value of the 
coupling g for metastability of currents with one unit of 
circulation, i.e., with 2nh/M , is (in the limit Rq 3> a osc ) 
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(8) 




Assuming the above factorization, one may start from the 
two-dimensional problem and derive the following one- 



i.e., it scales as w -1 / 2 . Clearly the numerical prefactor 
in the above expression depends on the form of the order 
parameter in the transverse direction R(p). If we define 

go = h 5 ^ 2 /(M 3 / 2 Wo' /2 i?o); this expression may be written 
in the more convenient form 



(9) 



Within this simplified model, the obtained power-law 
dependence of g m i n scaling as us^ 1 / 2 is in agreement 
with the numerical results plotted in Fig. 2. This figure 
shows that <7 m i n indeed increases linearly as a function of 
(uj/ujo)^ 1 / 2 for strong confinement in the transverse di- 
rection, when the motion is quasi-one-dimensional. From 
the numerical data we also see that for uj/ujo ^ 1, <? m in 
grows faster, as the extra degrees of freedom associated 
with the motion in the transverse direction start to play a 
role. We also notice that for small values of uj, when the 
confinement becomes weak, the maximum value of the 
confining potential at the center of the trap decreases. 
For example, for uj/ujo ~ 0.25, the density of the system 
is already substantial at the center of the trap, and the 
deviations from quasi-one-dimensional motion are very 
pronounced. 



B. Population imbalance and metastability 

We now turn to the question of metastability in a two- 
component system. Again, we study the minimum inter- 
action strength for the existence of metastable flow when 
the atoms are confined in a (two-dimensional) annular 
trap. To do this, we take Rq/ciq = 4 as before, and con- 
sider a fixed and relatively weak confinement uj/ujq = 1, 
varying xa = Na/(Na + Nb). This can be done eas- 
ily within the Gross-Pitaevskii scheme, where one can 
choose any value for the relative population Nb/Na, but 
not within the diagonalization approach, where one has 
to specify the (integer) particle numbers Na and Nb, and 
therefore the values of x a are more restricted. In addi- 
tion, in order to achieve small steps in xa, the number of 
particles has to be fairly large, which results in Hamilto- 
nian matrices of large size, making it hard to diagonalize 
numerically. 
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FIG. 2: (Color online) Minimum interaction strength <? m in/<7o 
for the existence of persistent currents, as a function of 
(ai/wo) _1//2 in an annular trap with Ro/ao — 4. The (black) 
circles show the result from the mean-field approximation, 
and the (red) crosses from the diagonalization of the Hamil- 
tonian. 
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FIG. 3: (Color online) The minimum interaction strength 
ffmin/ffo for stability of persistent currents in an annular trap 
with Ro/ao = 4 and u/uiq — 1, as a function of the relative 
population xa, obtained within the mean- field approxima- 
tion. The (black) circles correspond to L/N = 1, the (red) 
crosses to L/N = 2, and the (blue) squares to L/N = 3. 



As before, one can obtain a simple analytic result for 
<7min from the result of a strictly one-dimensional trap- 
ping potential Following the same arguments as in 
the previous subsection, then Eq. implies that, with 
the assumption of a step function for the density of the 
gas in the transverse direction, 
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which clearly reduces to Eq. (0) when xa = 1. According 
to this simplified model, g m - ln diverges at x a = 3/4 and 
thus metastability cannot exist for x A < 3/4. 

In Fig.[3]we show the numerical results that we obtain 
for the minimum interaction strength necessary for the 
existence of metastable states as a function of xa, within 
the mean-field approximation. The (black) circles corre- 
spond to a value of the angular momentum per particle 
L/N equal to unity. Our calculations indicate a diver- 
gence of g m in for sufficiently small values of xa, pretty 
much like the one-dimensional case, where this diver- 
gence occurs at xa = 3/4. On the other hand, metasta- 
bility takes place at a stronger interaction strength than 
in the purely one-dimensional case. Therefore, the de- 
viation from this limit works against the stability of the 
currents. This result is consistent with the general state- 
ment that metastability of the superflow is not possible 
in trapping potentials which decrease monotonically with 
the distance from the center of the trap [19j . 

What has been mentioned so far has to do with values 
of the circulation equal to one unit, i.e., equal to 2nh/M. 
According to Ref. 8], in the strictly one-dimensional 
problem persistent currents with circulation higher than 
one unit are stable only in the single-component case, i.e., 



for xa — 1 (or equivalently Nb = 0). Interestingly, un- 
like the one-dimensional case, our mean-field calculations 
show that persistent currents with a value of the circu- 
lation which is higher than one unit are stable for finite 
admixtures of a second component. This result is shown 
in Fig. 3 in the two higher curves, represented as (red) 
crosses, and (blue) squares, corresponding to L/N = 2 
and L/N = 3, respectively. 

The two limiting cases of strong and weak confine- 
ment are worth commenting on: as one approaches the 
one-dimensional limit, i.e., when u increases, the two 
curves corresponding to higher values of the circula- 
tion approach infinity more rapidly, eventually becoming 
vertical, as in the purely one-dimensional case. Thus, 
metastability with L/N — 2 or L/N = 3 is not possi- 
ble in this limit, in agreement with the results of Ref. 
Metastability is not possible, when the confinement be- 
comes very weak, either, i.e., when the width of the an- 
nulus gets large, and the gas has a significantly nonzero 
density at the center of the trap. In other words, the two 
curves have an infinite slope in the two limiting cases (of 
small and large values of uS), but have a finite slope for 
intermediate values, indicating the stability of persistent 
currents in this regime. 

While the behavior of the system looks similar in 
both limits, the underlying physical origin of the non- 
monotonic behavior of the slope of these curves is differ- 
ent in the two extremes. In the limit of small w, this is 
due to the fact that metastability is absent in an almost 
homogeneous gas [ll|. In the opposite limit of large to, it 
is the one-dimensionality of the problem that causes this 
effect @. 

The results that we have obtained from the numerical 
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FIG. 4: (Color online) The dispersion relation E(L) ob- 
tained from the diagonalization of the Hamiltonian, in a (two- 
dimensional) annular trap, with N = 17 atoms, uj/ujq = 1, 
Ro/ao = 4, and g/go = 30. The (red) crosses show the single- 
component case, with Na = 17 and Nb = (xa = 1), while 
the (black) squares correspond to a two-component system, 
with N A = 15 and N B = 2 (x A ~ 0.88). While in the first 
case there are two local minima in E(L) at L/N — 1 and 
L/N = 2, in the second case there is only one minimum, 
around L/N — 1. 

diagonalization of the Hamiltonian are shown in Fig. 0] 
and are consistent with those of the mean-field approxi- 
mation. Figure HI shows the dispersion relation, i.e., the 
lowest eigenvalue of the Hamiltonian for each value of L, 
with xa = 1 (red crosses) and xa ~ 0.88 (black squares). 
More specifically, we have considered N = 17 atoms, with 
Na — 17 and Nb = in the one case, and Na = 15 
and Nb — 2 atoms in the other one. In this calcula- 
tion the single-particle basis includes all the states with 
| to | < 5, and we consider the same parameters as before, 
with i?o/ao = 4, oj/loq = 1, and g/go — 30. 

The dispersion relation E(L) in the case of a sin- 
gle component has two local minima, at L/N = 1 and 
L/N = 2, showing that persistent currents are stable in 
both cases. This is consistent with the result obtained 
from the mean-field calculations shown in Fig. 3. On the 
other hand, in the case of two components, E(L) has 
only one local minimum around L/N = 1. Again, this 
is consistent with the result of Fig. 3, since for the corre- 
sponding values of xa and gmin/soj only the state with 
L/N = 1 is met ast able. 

V. SUMMARY AND CONCLUSIONS 

According to the present study, Bose-Einstein con- 
densed atoms which move in annular traps offer a very 
interesting system for the study of persistent currents. 

In principle one could investigate a surface in a three- 



dimensional phase diagram showing the minimum inter- 
action strength for stable currents as function of both the 
confinement strength u>, which determines the deviation 
from purely-one-dimensional motion, and the population 
imbalance xa ■ Given the size and the numerical effort of 
such a calculation, we have instead restricted ourselves to 
certain parts of this diagram, and in particular we have 
investigated two different questions. 

Firstly, we considered the case of one component, 
and showed that as the width of the annulus becomes 
larger, and thus there are deviations from the strictly 
one-dimensional motion, the critical coupling g m - ln that 
is necessary to achieve metastability increases. This re- 
sult is consistent with the fact that a necessary condition 
for metastability is that the confining potential - and as a 
result the particle density - does not decrease monotoni- 
cally from the center of the trap; in the picture of vortex 
dynamics, there is a force acting on the vortex that is 
associated with the rotational motion of the superfluid, 
which is in the opposite direction of the gradient of the 
single-particle density distribution. 

Secondly, we investigated the case of a mixture of two 
components, in a fixed annular potential. We found that 
ffmin which corresponds to persistent currents with one 
unit of circulation increases with the addition of the sec- 
ond component. Also, we found a similar behavior as 
in the one-dimensional problem, where g m i n diverges at 
xa = 3/4. Thus, the main effect of the finite width of 
the annulus is a relative increase in the critical coupling. 

On the other hand, a novel result that is absent in the 
one-dimensional case is the stability of persistent currents 
with circulation for the larger component which is higher 
than one unit. When the motion is one-dimensional, g m i n 
is infinite, and there are no stable currents. As the width 
of the annulus increases, these currents become stable for 
interaction strengths which are sufficiently strong, but 
finite. When the annulus becomes even wider, g m - ln be- 
comes infinite again. Therefore, the effect of the devia- 
tion from the one-dimensional motion is rather dramatic 
in this case. 

From the above results it is clear that the phase dia- 
gram g mm = g m i n (oj, xa) has an interesting structure and 
that the physics of persistent currents of Bose-Einstein 
condensed gases trapped in annular potentials is very 
rich. With the recent progress in building such trapping 
potentials in the laboratory, it will be of much interest 
to investigate all these effects also experimentally. 
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